Inelastic Scattering from Local Vibrational Modes 
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We study a nonuniversal contribution to the dephasing rate of conduction electrons due to local 
vibrational modes. The inelastic scattering rate is strongly influenced by multiphonon excitations, 
exhibiting oscillatory behaviour. For higher frequencies, it saturates to a finite, coupling dependent 
value. In the strong coupling limit, the phonon is almost completely softened, and the inelastic cross 
section reaches its maximal value. This represents a magnetic field insensitive contribution to the 
dephasing time in mesoscopic systems, in addition to magnetic impurities. 



PACS numbers: 73.23-b,73.63-b,72.10.Fk 

The loss of quantum coherence, occurring through in- 
clastic processes during which some excitations are left 
behind and the outgoing state is not the a single par- 
ticle state any more, is characterized by the dephasing 
time, tj,. The dephasing time can reliably be estimated 
from the low-field magnetoresistance data from the weak- 
localization corrections to the Drude conduct ivity[l|. In 
general, any dynamical impurity with internal degrees of 
freedom can dephase conduction electrons. If the impu- 
rity changes its state, the environment felt by the elec- 
trons changes, causing dephasing. 

Recently, dephasing due to magnetic impurities has 
been addressed successfully in the relevant temperature 
range (around and below the Kondo temperature Q) 
by relating the dephasing time to the inelastic cross 
section^, 0, @. The latter can be determined from 
the knowledge of the many-body T-matrix. It accounts 
nicely for various features of the dephasing time such 
as its saturation above the Kondo temperature. How- 
ever, in many systems, Coulomb-type interaction is not 
the only source of correlation, which leads us to con- 
sider the interaction of electrons with local vibrational 
(phononic) modes The interest to study these model 
is at least threefold: first, artificial quantum dots based 
on single moleculesQ (e.g. Cqo) usually distort upon 
the addition or removal of electrons, and can act as a 
quantized nano-mechanical oscillator by single electron 
charging^. Moreover, other realization contains sus- 
pended quantum dot cavity Q, where confined phonon 
modes influence the transport. Second, in many strongly 
correlated systems, such as heavy fermions or valence 
fluctuation systems, lattice vibrations are known to cou- 
ple strongly to electrons. The recent discovery of magnet- 
ically robust heavy fermion behaviour in filled skutteru- 
dite compound^ SmOs4Sbi2 renewed interest in Kondo 
phenomena with phononic origin 1^, U | • Both the spe- 
cific heat coefficient and coefficient of the quadratic tem- 
perature dependence of the electrical resistivity were 
found to be almost independent of an applied magnetic 
field. 

Finally, yet another motivation arises from recent ex- 
periments on Cug3Ge4Au3 thin films [12], where the ob- 



served dephasing time was found to be rather insensitive 
to high magnetic fields, suggesting non-magnetic dephas- 
ing processes, such as dynamical defects. 




FIG. 1: (Color online) Schematic picture of a nanomechanical 
resonator coupled to a conduction electron bath. 

Therefore, it is timely to investigate other source of 
dephasing, than magnetic impurities, such as a localized 
oscillator [6| . In general, these dynamic point defects are 
randomly distributed in a system. By assuming that the 
impurities are distributed sufficiently far from each other, 
it suffices to study one single mode, with fixed energy ujq 
(Einstein oscillator), similarly to the treatment of mag- 
netic impurities [2j. The phonon can couple to the local 
charge density, and modify the local potential felt by the 
electron (single impurity Holstein model (13. fli! llEl Tlijl). 
In the most general case, this amounts to study the Yu- 
Anderson or single impurity Holstein model (visualized 
in Fig. □} given by (Hill, 
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which describes 3-dimensional electrons interacting with 
a local bosonic mode, \&(r) denotes the bulk 3D elec- 
tron field operator and Cfc its Fourier components. This 
accounts for intrinsic and extrinsic point defects such 
as interstitial or substitutional impurities (2fj|. From the 
latter class, foreign atoms (like Fe in Ag or Au[3, H|) 
mainly represent a source of potential scattering. Be- 
yond this, they are capable of magnetic scattering 0, 0], 
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and can cause a modification of the local phononic en- 
vironment as well, resulting in a Hamiltonian like Eq. 
(fl]). Besides foreign atoms, the self- interstitial impurities 
(proper atom at a non-lattice site) also distort the local 
potential background felt by the conduction electrons, 
leading to local electron-phonon interaction. This model 
has a rich history in the field of nonlinear impurities 21 1 . 
In quantum optics, such model arises in arrays of non- 
linear optical waveguides for example. Moreover, in the 
condensed matter context, Eq. (TTJ) bears close resem- 



blance to the Anderson-Holstein impurity model[22j, |23j, 
describing molecular vibrations in quantum dots. It can 
be regarded as a generalization of scattering on two-level 
systems to the multi- level case[24j]. 

Due to the isotropic nature of the electron-phonon cou- 
pling, the electron field operators can be expanded in the 
appropriate angular momentum channels (depending on 
the dimensionality of the conduction electrons) , and only 
the s-wave component couples to the impurity [25| as 
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where R(x) stands for the chiral, right moving radial 
component, r > is the radial coordinate, hp is the 
Fermi wavenumber. Thus, the model can be mapped 
onto an effective model of chiral fermions, interacting 
with the bosonic mode at the origin, and is ideal for 
bosonization (R(x) = exp [i-\/47r$ {x)]/V 2wa) : 



H = - 



oo 

f P 2 mw 2 

iv / dxR+(x)d x R{x) + gQp(0) + — + — ^Q 2 , 
J 2m 2 



(3) 

and only the radial motion of the particles is accounted 
for by chiral (right moving) fermion field, p(x) —: 
R + (x)R(x) : is the normal ordered electron charge den- 
sity, v is the Fermi velocity, g — gdkp/ir describes 
the local electron-phonon coupling, m and ojq are the 
phononic mass and frequency, respectively, Q and P are 
the phonon displacement field and momentum conjugate 
to it. In the spirit of the Yu- Anderson model[18(, the 
phonons couple to the electronic density fluctuations in 
our model as well. The ground state will be degenerate, 
similarly to the Kondo problem between occupied and 
unoccupied electron states, leading to the formation of a 
double well potential[l8|. The bosonized Hamiltonian is 
identical to that of coupled harmonic oscillators 2(| 27 1 . 

Due to the electron phonon coupling g 1 the phonon 
mode softens as|17l \21\ 



cj p± = -»r ± ^uo 2 - r 2 - r^ 2 /r 2 , (4) 

where T 2 = ttoj 2 /4W <C ujq <C W, W is the bandwidth of 
the conduction electrons, and T — n(gp) 2 /2m for small 
g, and approaches r 2 as g — > oo. Here, p — 1/2ttvf is 



the chiral electron density of states. The explicit depen- 
dence of r on g cannot be determined by the bosonization 
approach [2?]] ■ The real part of the phonon frequency re- 
mains finite (underdamped) for T < T% w r 2 (l — T 2 /lli 2 ). 
For Ti < T < r 2 , the oscillatory behaviour disappears 
from the phononic response (~Reu> p ± — 0), and two dis- 
tinct dampings characterize it (overdamped) . This soft- 
ening is observable in many physical quantities, such as 
the enhancement of the local charge response, indicating 
the charge-Kondo phenomenon. 

The local retarded Green's function of the conduction 
electrons is defined, using Eq. ([2]) aspjj 

G R (t) = -ie(t)({^(x = 0,t),^ + (y = 0,t = 0)}) = 
= -m)f- E lim({i?(7a,i),i? + (7'«,0)}) = 

47T ^ — 4 Q— >0 

7,7'=± 

= -'9W§ E (exp[C 7 ,y(*)]+exp[CV, 7 (-<)]), 



7,7'=± 



(5) 



where C 7n >(t) = lim Q ^ ^(^(-ya, t)$(-f'a, 0) - 
[$(70;, t) 2 — $(7'a, 0) 2 ]/2) — ln(a), whose second appear- 
ance in Eq. (O follows from time reversal symmetry, 
and a ~ 1/W is the short distance cutoff in the 
bosonized theory, 7 and 7' denotes the sign of the space 
coordinate. The impurity site is defined in terms of the 
chiral operators as the limiting value from Eq. ([2|). The 
7 = 7' combinations does not contain any information 
about the presence of the impurity at the origin[29j]. 
since both field operators are taken on the same side 
of the impurity, and no scattering is experienced. 
On the other hand, (R^a, t)R + (— 7a, 0)) describes 
scattering through the impurity, and depends strongly 
on its properties. The expectation value C 7i7 '(i) at 
bosonic Matsubara frequencies can be evaluated from 
the bosonized Hamiltonian [27| . Then by analytical 
continuation to real frequencies, we can determine the 
spectral intensity following Ref. H^, which leads to the 
desired function, C 7i7 '(i). This follows the derivation 
of the position autocorrelator of a harmonic oscillator 
coupled to a heat bath[26j]. The local retarded Green's 
function follows as 



G r (lo) = -tir^ 1 + F(«) 



4Ti/(u> p+ — cj p _) n 



(6) 



where 



F(w) 



exp(iwi) 
at lm exp 

7rf 



4r(c?[iii> p+ ] — d[tuj p -] 



i{tjj p+ - uj p ^) 



(7) 



and d[x] — exp(-ix) (— Ci(— x) + i (Si(— x) — f )) in the 
underdamped case, Si(x) and Ci(a;) are the sine and co- 
sine integrals, and pb — k 2 F /2-K 2 VF is the 3D bulk density 
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of states. When Eq. © is expanded in T, it agrees 
with perturbative results[31|. We can identify the T- 
matrix using the Dyson equation for the local Green's 
function as G(0,w) = G o (0,w) + G (0,w)T(w)G Q (0,u), 
where Go(0, uj) = —iirpi, is the unperturbed Green's func- 
tion. Therefore, we determine the full T-matrix as 



1 u)/u)o 2 



30 



3» = 



2irp b 



4ri/(w p+ -cjp_) 



+ F(w)-1 (8) 



which is the main result of this paper. It contains all the 
non-perturbative effects of the quantum impurity intro- 
duced in Eq. ((TJ) . As was shown in Ref . d, this possesses 
all the information needed to evaluate the various cross 
sections. From this, we obtain the many-body s-matrix 




FIG. 2: (Color online) The evolution of the s-matrix on the 
complex plane is shown at T = for W = Wljo, r/r 2 = 
0.3, 0.5, 0.7 (blue dashed), 0.8 (blue lines), 0.9 (green line), 
0.94 (red dashed line) and 0.98 (red solid line) from right to 
left. Close to the critical coupling F = r 2 , the scattering is 
dominantly elastic for small frequencies, and follows the unit 
circle. Then, it reaches rapidly the maximal inelastic rate 
(s = 0) with increasing w, because at the critical coupling, 
the phonon excitation energy \u)p-\ — > 0, facilitating inelastic 
processes. 



as[32[ 



s(u>) = 1 — i2ir pi,T (u>) , 



(9) 



which stays always within the unit circle in the complex 
plane. The circumference of the unit circle denotes per- 
fect elastic scattering. The total, elastic and inelastic 
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FIG. 3: (Color online) The frequency dependent scattering 
cross sections are shown for W = lOwo, r/r 2 = 0.3, 0.6, 0.9 
and 0.99 from left to right, top to bottom. The solid blue, 
red dashed and black dashed-dotted line denotes a tot, <Jinei 
and CT e i, respectively. For small frequencies, all scattering 
rates vanish. Then, wild oscillations show up, representing 
the various multiphonon (n-phonon) processes, setting in at 
n x At higher frequencies (|Wp±| <C w <C W), in- 

elastic processes dominate. With increasing Y, the weight of 
both elastic and inelastic processes increase. By approaching 
the critical coupling, the inelastic scattering rate reaches its 
maximal value, ctq /4, because with the vanishing phonon fre- 
quency, uj p -, phonon excitations become more likely, favour- 
ing inelastic scatterings. Note the different horizontal and 
vertical scales on the figures. 



cross section is readily evaluated from s(lu) as0,0, HH 

(10) 

(11) 
(12) 



<Jtot = y-(i ~ R es M)- 

CineZ — 
a , 



l*Hl a ), 



a el = -^(l-2Re»(a;) + |s(w)n, 

where (To = 2/7tp u v is the maximal value of the cross sec- 
tion. (Tinei is directly related to the dephasing time t$, 
measured in weak- localization experiments 3] . A local- 
ized vibrational mode contributes to the inelastic scat- 
tering rate even at T = 0. The evolution of its s-matrix 
is shown in Fig. [2] In the weak coupling limit (small 
r), scattering is negligible. For small frequencies, the 
s-matrix follows the upper semicircle, thus scattering is 
mostly elastic. For larger T, the weight of the inelas- 
tic component increases with T until it reaches its max- 
imum value at s = at the critical coupling T2. There, 
the total scattering cross section is divided equally be- 
tween inelastic and elastic processes with weight cfq/A 
for \uj p ±\ <C UJ <C W, similarly to the channel isotropic 
two-channel Kondo model at ui = 0, which is a prototype 
model for non-Fermi liquid behaviour. However, the clas- 
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sification of impurity states as Fermi or non-Fermi liquid 
becomes inappropriate in the present case, since scatter- 
ing is caused by bosonic, and not fermionic modes. Since 
the transition occurs at g — > oo (r — * r 2 ), the crossing 
of this critical point is impossible. 

Fig. [3] shows the total, elastic and inelastic scat- 
tering rates as a function of frequency for several T. 
At zero frequency, all rates vanish, since the electron 
Green's function is pinned to its non-interacting value 
due to some generalized Fermi liquid relations [3, EH- 
For small couplings and u> <Rew p ±, the ground state 
contains no phonons but a filled Fermi sea, thus elec- 
trons scatter off elastically. By exceeding this threshold 
(Reu> p ±), the total scattering is mainly given by the in- 
elastic cross section, which develops steps at at n x \u> p ±\ 
(with n integer), stemming from n-phonon excitations. 
At n = 1, this is also predicted by perturbation theory [3l| 
as Otot ~ &inei = <JoTitQ(uj — u>o)/u>o- Higher steps occur 
with smaller weight in alternating fashion. At the same 
time, the elastic rate exhibits a small cusp at loq due to 
the logarithmic singularity of the real part of electron self 
energy (3lj|. 

For larger couplings, these sharp steps become smeared 
and turn into oscillations due to the finite lifetime of the 
phonons. These oscillation are not of Friedel type, cor- 
responding to the interference of incoming and scattered 
waves, but arise from multiphonon excitations. The elas- 
tic scattering rate grows progressively, and becomes the 
dominant scattering mechanism for 1 and 2 phonon ex- 
citations. 

By approaching the critical coupling T2, the inelas- 
tic scattering reaches its maximum, since the vanishing 
phonon energy, cj p _, makes arbitrary multi-phonon pro- 
cesses possible. We estimated, H| luq to fall into the 0.5- 
50 K range, which compares well to the typical tempera- 
tures of weak localization experiments. These results at 
T = are in marked contrast to those found for Kondo 
impurities d, 0, [H . 

In actual experiments (on e.g. Ag or Au samples with 
Fe impurities), the dephasing time is measured at fi- 
nite temperatures, and can be decomposed as l/r^ = 



l/ie-e + l/ r e— ph + l/ r imp- I n the high temperature 
limit, it is dominated by the bulk electron-electron 
(l/r e _ e ~ T 2 / 3 ) and electron-phonon (l/r e _ pft ~ T 3 ) 
interactions, and is accounted for by the Altshuler- 
Aronov-Khmelnitzky theory [3^|. At lower temperatures, 
the impurity contribution becomes important (l/ri mp ~ 
n-impVinei) , and scales with the impurity concentration 
riimp. Magnetic impurities account successfully for this 
latest contribution down to O.ITk (with Tjc the Kondo 
temperature). 

Our calculation on the local oscillator model can be ex- 
tended to finite temperatures. Conformal invariance can 
only be used for times much bigger (25| than 1/wq, which 
translates to T C wo, leading to T 2 corrections. The 
same conclusion borns out from the generalized Fermi 



liquid relations^, 14|. In the 4Reuj p ± < T <C W range, 
the inelastic scattering rate reaches its maximal value, 
Cinei = Oei = o"o/4, since at high temperatures, the 
phonon state is highly populated, facilitating inelastic 
scattering. This behaviour is similar to the T = 
case close to the softening r — > T2, where the vanish- 
ing phonon excitation frequency favours inelastic pro- 
cesses. In between these two regions (T ~ too), a smooth 
crossover takes place. 

The Kondo picture starts to deviate from the exper- 
imental data below the Kondo temperature 0, [U l35| . 
This mismatch is cured by introducing a very small 
amount of magnetic impurities with a very low Kondo 
temperature [B[ or by adding a constant background [35| . 
Instead, a flat, temperature independent dephasing time 
(with ujo -C Tk <C W) follows naturally from local vi- 
brational modes, as was demonstrated above, and serves 
as an obvious explanation for the measured behaviour. 
Therefore, local vibrational modes can account succes- 
sively for the measured dephasing time in weak localiza- 
tion experiments. 

The second crucial difference with respect to Kondo 
impurities arises from the magnetic field insensitivity of 
the phonons, and requires further experiments to reveal 
the differences [l2j|. 

In conclusion, we have demonstrated through an exact 
solution of the single impurity Holstein or Yu- Anderson 
model, that local vibrational modes can have a strong 
impact on the dephasing time of electrons. The inelas- 
tic scattering rate exhibits strong oscillations at frequen- 
cies comparable to the phonon excitation energy, and 
then saturates to a finite, coupling dependent value. At 
the extreme strong coupling limit, close to the complete 
softening of the phonons, the s-matrix vanishes and the 
inelastic cross section reaches its maximal value. This 
phonon mediated scattering mechanism is expected to be 
rather insensitive to the applied magnetic field, in con- 
trast to Kondo-type impurities, and can contribute to 
the dephasingtime in certain alloys containing dynami- 
cal defects 
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